Domain wall motion in ferromagnetic nanowires driven by arbitrary time-dependent 

fields: An exact result 
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We address the dynamics of magnetic domain walls in ferromagnetic nanowires under the influence 
of external time-dependent magnetic fields. We report a new exact spatiotemporal solution of 
the Landau-Lifshitz-Gilbert equation for the case of soft ferromagnetic wires and nanostructures 
with uniaxial anisotropy. The solution holds for applied fields with arbitrary strength and time 
dependence. We further extend this solution to applied fields slowly varying in space and to multiple 
domain walls. 
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Introduction. — The motion of magnetic domain walls 
(DWs) in ferromagnetic nanowires has recently become 
a subject of intensive research in the condensed matter 
physics community [l[ . Manipulation of DWs by external 
magnetic fields, and in particular, the question of how the 
DW propagation velocity depends on the applied field 
have drawn considerable attention 0-0] • 

In ferromagnetic nanowires, the dynamics of the ori- 
entation of the magnetization distribution, m(x, t) (nor- 
malized so that |m| = 1), is described by the Landau- 
Lifshitz-Gilbert (LLG) equation 5] 



9m dm , 

— +«mx — = (1 



(H(m) + H Q 



(1) 



where x is the coordinate along the nanowire, t is time, 
a is the Gilbert damping parameter, H a denotes the ap- 
plied magnetic field, and H(m) = —SE/Sm, where 



^ = 4/ E ^ ^ + f/ R (l- (m .*f)d, 



K2 , , 
-7T / (m-y) 



Ax. 



(2) 



is the reduced micromagnetic energy. Here, A is the ex- 
change constant of the material, and K\ , K 2 > are the 
anisotropy constants along the (easy) x- and (hard) y- 
axes. The anisotropy constant along the z-axis is taken 
to be zero by convention. 

To date only one exact spatiotemporal Q solution of 
the LLG equation has been reported in the literature, 
namely the so-called Walker solution [6] . The analysis of 
Schryer and Walker [6] applies to the case where K 2 > 0, 
ie where the anisotropy constants in the transverse plane 
are strictly unequal. This is appropriate for a thin film 
or thin strip geometry. The applied field is taken to be 
uniform in space, constant in time, and directed along 
the nanowire, i.e., H a (x, t) = -ff a x. For |if | less than 
a certain threshold Hy/ , the so-called Walker breakdown 
field, a planar domain wall propagates rigidly along the 
nanostrip with velocity depending nonlinear ly on H a . 

In this Letter we present an exact spatiotemporal so- 
lution of the LLG equation that, to our knowledge, has 



not been previously reported in the literature. We con- 
sider the case of transverse isotropy, ie K 2 = 0. This is 
appropriate for soft ferromagnetic nanowires whose cross- 
sectional dimensions are comparable, as well as for uni- 
axial nanowires whose easy axis lies along the wire. We 
take the applied field to lie along the nanowire, as in the 
case of the Walker solution, but allow for arbitrary time 
dependence, i.e., H a (x,t) — H a (t)x. 

Exact solution of the LLG equation. — The boundary 
conditions appropriate for a domain wall with finite mi- 
cromagnetic energy E(m) are given by m(x, t) —y ±x as 
x — > ±00. For K 2 = the magnetization-dependent field 
H is given by 



H(m) = A 



d 2 m 

dx 2 



K±(m ■ x)x . 



(3) 



We now take into account the fact that m has its values 
on S 2 , and parametrize m in terms of angles 0(x,t) and 
4>{x,t) according to m = (cos 9, sin 9 cos <f>, sin 9 sin <fi). 
From Eqs. (JTJ) and © we obtain the LLG equation in 
the equivalent form 



9 - a<j> sin 9 + A(l + a 2 ) {<t>" sin 9 + 26' (p' cos 9) = , 



(4a) 

a9 + 4> sin 9 + (1 + a 2 ) ( - AO" + A(cj)') 2 sin 9 cos 9 

+ Ki cos 9 sin 9 + H a (t) sin 9) = , (4b) 

where dot ' denotes d/dt and prime ' denotes d/dx. 
We now look for a solution of Eq. ^ in the form 



,(x,t) = 9 (x - £*(*)) , 4>*(x,t) = , (5) 



where 



%{x) = 2 arctanexp (— x/do) , do = yA/Ki 



(6) 



9q(x) describes the static domain wall in the absence of 
an applied field. The magnetization density determined 
by 9o(x) and <fio(z) = tt/2 minimizes the micromagnetic 
energy E(m) for the specified boundary conditions. Sub- 
stituting Eq. ([6]) into Eq. (j4]), and taking into account 



2 



that 6q = — sin 8 /do and 0' ' 
that 9„ and 
that x*(£) and 



sin26» /(2^), we find 
satisfy the LLG equation (|4| provided 
satisfy 



= -ad H a (t) , 



-#«(*)• 



(7) 



(In fact, © and (0 provide the only solution of the form 
©•) 

Equations ([5][7]) constitute the main result of this Let- 
ter. They represent an exact solution of the LLG equa- 
tion, and describe a DW, with profile independent of the 
applied field, propagating along the nanowire with veloc- 
ity ±* while precessing about the nanowire with angular 
velocity No restrictions have been imposed on the 
strength of the applied magnetic field and no assumptions 
have been made about its time dependence. 

We now compare the precessing solution Eqs. ([5][7]) 
with the Walker solution 6]. The Walker solution is de- 
fined only for K 2 > (the fully anisotropic case) and 
time-independent H a less than the breakdown field 



H w = aK 2 /2. 



It is given by 



9 w {x,t) = 6 ( - — I , <p w (x,t) 



(8) 



<t>w , (9) 



where 



sm.2(j> w = H a /Hi 



(10) 



determines the (fixed) inclination of the DW plane and 



1 + a 2 

Vw = 7 d H a , 7 = 



K\ + K 2 cos 2 cf>w 



1/2 



(11) 



gives the DW velocity. 

There are several characteristic differences between the 
Walker solution and the precessing solution which should 
be distinguishable experimentally. Foremost is the fact 
that the Walker solution exists only for constant applied 
fields whose strength does not exceed a certain thresh- 
old, so that the DW velocity is bounded. The precessing 
solution is defined for time-dependent applied fields of ar- 
bitrary strength, so that the DW velocity, which for the 
precessing solution is proportional to the field strength, 
can be arbitrarily large. Next, while for the Walker solu- 
tion the plane of the DW remains fixed, for the precessing 
solution it rotates about the nanowire at a rate propor- 
tional to H a . Finally, we observe that, for the Walker 
solution, the DW profile contracts (7 > 1) or expands 
(7 > 1) in response to the applied field, whereas for the 
precessing solution the DW profile propagates without 
distortion. 

Spatially nonuniform applied fields and multiple do- 
main walls. — We now extend our results to applied fields 
that depend on both position along the nanowire and 



time, i.e, H a = H a (x,t)x. For any (non-singular) ap- 
plied field, Eq. (|4]) is satisfied at x outside the DW tran- 
sition layer \x — x*(t)\ do (up to exponentially small 
terms). Assuming now that the field varies slowly across 
the transition region, 

\H a (x,t) - H a (x*(t),t)\ < \H a (x*(t),t)\ 

for \x-x*(t)\ < d , (12) 

we obtain an approximate solution of the LLG equation: 
the magnetization density is given by Eqs. ((5|) and ([6|) 
with 

x* = — adoH a (x M (t),t) , 0* = -H a (x*(t),t) . (13) 

The physical meaning of Eq. (113[) is quite obvious: the 
DW is only sensitive to the applied field within the tran- 
sition layer. 




FIG. 1: (Color online) Dynamics of domain walls. See text 
for discussion. 

This approximation can now be extended to the case 
of N non-overlapping DWs. Indeed, 



N 



e N (x,t) = ^0 o ((-i) n+1 (z - x n {t))) 



(14a) 



4>N{x,t) = 4>n(t) , n = n minimizes \x — x n (t)\ , (14b) 

with Xk+i(t) — Xk{t) 3> do for k = 1, . . ., N — 1, consti- 
tutes an approximate solution of the LLG equation given 
that 



x n = (-l) n ad H a (x n (t) 1 t) 
0n = -H a (x n (t),t) , 



(15a) 
(15b) 



for n = 1, ...,N. For the case of a spatially uniform 
applied field Eqs. (TT41 and (|T5|) describe the time evolu- 
tion of N DWs such that any two adjacent WDs travel 
in opposite directions while rotating in the same direc- 
tion (and with the same angular velocity) around the 
nanowire. 

Conclusions. — In this Letter we have presented an ex- 
act spatiotemporal solution of the LLG equation that has 
not been previously reported in the literature. The va- 
lidity of the new solution requires no assumptions about 
the time-dependence or strength of the applied field. 

We have then provided a natural extension of the so- 
lution to physical situations in which the applied field 
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varies slowly in space. An approximate solution of the 
LLG equation for the case of multiple domain walls has 
also been obtained. 
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